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INTRODUCTION 


1 . Collaborators in Japan for MAGSAT Data Analysis 

The Japanese MAGSAT Team (the initial membership is shown on the cover 
page) carried out investigations with MAGSAT data over the past 28 months 
with the assistance of the following collaborators. 

I. Nakagawa (Earthquake Research Institute, University of Tokyo, Tokyo 113) 
T. Nakatsuka and Y. Ono (Geological Survey of Japan, Tsukuba 305) 

M. Yanagisawa (Institute of Space and Astvcnautical Science, Tokyo 153) 

T. Kamei, T. Iyemori, S. Tsunomura and T. Kumaki (Geophysical Institute, 
Kyoto University, Kyoto 606) 

II. Ejiri and H. Sakurai (National Institute of Polar Research, Tokyo 173) 

The following persons participated in the production or analysis of ground 
or airborne magnetic surveys, the data of whicn will be fully used in the 
near future for comparison with MAGSAT data. 

Y. Ueda, K. Onodera, T. Kaneko and R. Nakagawa (Hydrographic Department, 
Maritime Safety Agency, Tokyo ^04) 

S. Fujita (Magnetic Observatory, Kakioka, Ibaraki-Pref . 315-01) 

K. Shibuya (National Institute of Polar Research, Tokyo 173) 


2 . Advisory Committee to the MAGSAT Team 

In September 1979 an advisory committee for the MAGSAT Team was set up 
by the National Committee for IUGG, with the purpose of (i) giving appropriate 
scientific advice to the MAGSAT Team, and (ii) promoting the future possible 
utilization of MAGSAT data in the various disciplines beyond the subjects 
chosen by the MAGSAT Team for the immediate study. The MAGSAT Team welcomed 
the participation of the advisory members in the meetings held every four 
months throughout the investigation period. 


3. Financial Support to the MAGSAT Team 

In addition to financial support from each participating institution, 
the team could receive a special grant-in-aid from the Ministry of Education, 
Science and Culture for 1980-82, which was mainly used for (i) the meetings 
of the team held every four months, (ii) preparation and printing of the 
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Progress Reports, and (iii) reproduction and/or editing of MAGSAT data tapes 
for the workers in various locations in Japan. The Ministry, which endorsed 
the NASA-Japan cooperation of MAGSAT investigations, provided through the 
Institute of Space and Astronaut ical Science the expenses for participating 
in the MAGSAT Investigators' Meetings convened by NASA. 


A. Rese a rch Memoranda from the Japanese MAGSAT Team Colleagues 

In response to the request from the chairman of the Japanese MAGSAT 
Team, che team colleagues presented their research memoranda listed below. 


Memo No. 


Authors and Titles 


Al-1 Nakagawa, I. and T. Yukutake, Spatial properties of the geomagnetic 
field in the area surrounding Japan. 

Al-2 Ueda, Y. , T. Kaneko, H. Nakagawa, K. Onodera, and S. Oshima, 

Evaluation of Magsat field model MGST 4/81 and the regional 
geomagnetic anomaly around Japan estimated from airborne 
magnetic survey. 

A2-1 Nakatsuka, T. and Y. Ono, Geomagnetic anomalies around Japan 
derived from MAGSAT data. 

A2-2 Nakagawa, I. and T. Yukutake, Local geomagnetic anomalies derived 
from MAGSAT data over the area of Japan and its surrounding 
seas by use of Fourier analysis technique. 

A2-3 Fujita, S., M. Kawamura, T. Yukutake, Y. Ueda, T. Kaneko and 
S. Oshima, Regional magnetic anomaly around Japan. 

A3-1 Kaminuma, K. and K. Shibuya, Aeromagnetic survey in Antarctica. 

Bl-1 Yanagisawa, M. and N. Fukushima, External field correction for 
satellite magnetic data. 

Bl-2 Yanagisawa, M. and N. Fukushima, Dawn/dusk asymmetry of the ring 
current field observed by MAGSAT. 

Bl-3 Yanagisawa, M. and N. Fukushima, Eastward ring current inferred 
from MAGSAT and 0G0 5 data. 

Bl-4 Suzuki, A., T. Kamei, T. Kumaki and N. Fukushima, Electric current 
in space below the MAGSAT level. 

Bl-5 Maeda, H., T. Kamei, T. Iyemori and T. Araki , Geomagnetic perturb- 
ations at low latitudes observed by MAGSAT. 

Bl-6 Araki, T. , T. Iyemori, S. Tsunomura and T. Kamei, Sudden commence- 
ments of magnetic storms observed by MAGSAT. 

B2-1 Iijima, T., T. Araki and N. Fukushima, Large-scale field-aligned 
currents over the polar cap with reversed polarities and S^ 
associated with northward interplanetary magnetic field. 


The outlines of all these memoranda are included in this Draft Final Report, 
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and they are cited with their Memo. No. They were first checked and summa- 
rized by T. Yukutake (for solid-earth subjects) and T. Iijima (for external 
magnetic field subjects) before the present Draft Final Report was written 
by the chairman. 

We hope to publish most of the above memoranda prepared by the Japanese 
MAGSAT collaborators in early 1983, possibly in bulk in the "Journal of Geo- 
magnetism and Geoelectricity", as soon as such a publication is agreed upon 
by NASA. 

5. Additional Remarks 

The following 60 pages (27 pages of text, and 33 pages of 
figures) is a summary of the MAGSAT investigations carried out in Japan till 
November 1982. In the special issue for "MAGSAT Preliminary Results (Geo- 
physical Research Letters, Volume 9, Number 4, April 1982) the following five 
papers appeared as contributions from the Japanese MAGSAT Team, i.e. 

Yanagisawa, M., M. Kono, T. Yukutake and N. Fukushima, Preliminary inter- 
pretation of magnetic anomalies over Japan and its surrounding area, 
pp. 322-324 (Paper 2L0022). 

Maeda, H., T. Ivemori and T. Araki, New evidence of a meridional current 
system in the equatorial ionosphere, pp. 337-340 (Paper 2L0121) . 

Araki, T., T. Ivemori, S. Tsunomura, T. Kamei and H. Maeda, Detection of 
an ionospheric current for the preliminary impulse of the geomagnetic 
sudden commencement, pp. 341-344 (Paper 2L0157). 

Suzuki, A. and N. Fukushima, Sunward or anti-sunward electric current in 
space below the MAGSAT level, pp. 345-347 (Paper 2L0116). 

Iijima, T., N. Fukushima and R. Fujii, Transverse and parallel geomagnetic 
perturbations over the polar regions observed by MAGSAT, pp. 369-372 
(Paper 2L0355) . 

The contents of these published papers are not fully reproduced in this Draft 
Final Report, however the important conclusions in these papers are briefly 
included in connection with more recent results of analyses. 

In addition to the work in Japan, Y. Kamide has been cooperating with 
some U.S. scientists (during his stay in NOAA, Boulder, Colorado) in the 
comparison of MAGSAT data with other relevant uata (such as electron data by 
TIROS-N and NOAA-6). His results are not included in this report, but they 
will be known later through some papers by the investigators concerned in- 
cluding Y. Kamide. 
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SUMMARY OF THE RESULTS OF MAGSAT INVESTIGATION IN JAPAN 


A. CRUSTAL STRUCTURE NEAR JAPAN AND ITS ANTARCTIC STATIO N 

The objective is to study the crustal structure near Japan and its 
Antarctic base by constructing a model of the regional magnetic field and 
investigating the local magnetic anomalies and their origin 


A-l . Regional Magnetic Charts 

It has been attempted to construct a model of the regional geomagnetic 
field over the Japanese islands and their neighborhood for the epoch 1980.0, 
by means of MAGSAT data (Nakagawa and Yukutake, Memo. Al-1) and the 
aeromagnetic survey data obtained by the Hydrographic Department in 1980 
(Ueda et al., Memo. Al-2) . It was originally planned to make a model by 
incorporating the MAGSAT data with the surface data, with downward continu- 
ation of the MAGSAT data. incorporation of the MACSAT data is, however, 
left for future work. The field synthesized by the model of the third 
degree polynomials was compared with the existing reference field, IGRF 
1980.0 and MGST(4/81). Better agreement is seen with MGST(4/81) rather 
than with IGRF 1980.0 (Ueda et al.. Memo. Al-2). This seems to arise from 
the difference in the number of terms included in the models. Poor appro- 
ximation of IGRF is perhaps due to lack of terms higher than degree 10 
which also represents contribution of the cort field. 


1 . Spatial properties of the regional geomagnetic field 

For the magnetically quiet period of K^J-2^, three component data 
X, Y and 7 . from MAGSAT Investigator-B tapes along individual paths over 
the area of 120°-160°E and 18°-58°N were subjected to the polynomial fit, 
by means of Legendre functions, by varying thp maximum degree of series N 
from 1 to 9 (Nakagawa and Yukutake, Memo. Al-1). The study area was divided 
into subdivisions, 5 regions for the data of the ascending path and 5 for 
those of the descending path, in a strip bounded by orbital paths of every 
10° intervals of longitude. Fig. Al-1 shows the decrease in the magnitude 
of the coefficient with the degree of the polynomial terms for the Z 
component along descending paths in region 8, which covers the area from 
the southwestern part of Japan to the east coast of Siberia. The magnitude 
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of the coefficient decreases rapidly from degree 1 to 6, and beyond degree 
6 the rate of decrease becomes small. Similar features are seen with X 
and Y components and those in all other regions. If the magnitude of the 
coefficient is regarded as repres< nting the amplitude spectrum, this 
indicates that there are two distinct structures in the spatial spectrum 
of the observed field. One is the lower degree terms of large amplitude, 
and the other is the higher degree terms of small amplitude with low 
decreasing rate. 

The difference in the spectrum between the higher degree terms and 
the lower ones is more clearly seen in the residual field after removing 
the fitted polynomials. An example in Fig. Al-2 shows the rms residuals 
of the X-component plotted against the maximum degree of fitness for 

region 3, which covers the region from southwestern Japan to northeastern 
China through Korea. The rms residuals decrease sharply from degree 1 
to 5, and become almost constant beyond degree 5. Similar feature are 
observable with the Y and Z components as well. 

The above feature in the spatial spectrum is very similar to the 
spectrum structure obtained from the spherical harmonic analysis of the 
global data, suggesting that these structures are of the origin of the 
core and the crust respectively. The two types of field become comparable 
at degree 13, corresponding to the wavelength of about 3000 km (Langel 
and Estes, 1982). 

In the present analysis, however, the overlap of the spectrum of the 
core and the crustal field does not always occur at the same degree. The 
overlap degree for the individual orbit paths was averaged in each region, 
and the frequency distribution of the average overlap degree in the 10 re- 
gions is shown in Fig. Al-3. For the X component, the overlap degree is 
either 5 or 6 in 9 regions. For the Y component the critical degree 
concentrates at 3. On the other hand, for the Z component it spreads 
toward higher harmonics up to degree 7. 

From the degree of the polynomial one may roughly estimate the wave- 
length of the field. The degree n was related to wavelength A by the 
equation 

\ = L/[(n-l)/2] 

where L is the length of the path used for the analysis, which is taken 
as 7000 km. The overlap degrees of 5 and 6 correspond to the wavelength 
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of 3500 to 2800 km. This is in agreement with the overlap wavelength 
3000 km derived from the spherical harmonic analysis of worldwide data. 

2 . Regional magnetic chart s 

Aeromagnetic purveys were conducted by the Hydrographic Department of 
the Maritime Safety Agency in 1980 with fluxgate magnetometers for the 
three components (Hydrographic Department, 1982). Accuracy of the measure- 
ment was about 10 nT for the Z component, and somewhat worse for the other 
components, mainly controlled by the accurac} of orientation <-f the instru- 
ment. Total intensity was computed from the component data. Fig. Al-4 
shows the flight tracks; the mean altitude of the flight was 3200 m. The 
surveyed data along the flight tracks were subjected to the running mean 
over 60 km intervals, and expanded in a polynomial series of the third 
degree with longitude and latitude, as 

F, X, V, Z = a l + a 0 $' + a^X' + a^' 2 + a^'\' + a^’ 2 
+ a^' 3 + a^' 2 \' + a g *'\'~ + a 1Q \' 3 

where ;'(in radian) = ^(latitude)-^, and X'(in radian) = \ (longitude)- \ Q , 

with .? = 34.0°N and \ =136.0°E. The coefficients a., a„, a._ were 

o o 1 1 ID 

determined by a least squares method so as to satisf\ the conditions 

-* 2 2 2 2 

simultaneously (rot F ) ^=0 , and F =X +Y~+Z (Ueda et al.. Memo. Al-2) . The 
calculated coefficients are listed below. 



1 F 1 

X 

Y 

Z 

a l 

45992.93 

31272.59 

-3376.88 

33555.50 

a 2 

31693.25 

-24319.86 

-7810.60 

65319.69 

a 3 

-16232.68 

-4586.95 

3026.70 

-17669.88 

a 4 

8583.76 

-17548.71 

8402.03 

-29365.45 

a. 

a 

-5192.29 

25466.02 

-6962.53 

-15106.24 

a, 

0 

-3335.83 

-3732.36 

15062.43 

-754.07 

a, 

/ 

-43638.32 

28637.06 

-7738.63 

-32772. S6 

a 8 

46682 . 60 

-24573.14 

5665.64 

106022.66 

a 9 

-20867 . 56 

7335.80 

-11165.48 

-38778.18 

a 10 

23507.17 

-5893.14 

0.0 

37060.18 
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Figs. Al-5(a) to (d) show the computed values of F, X, Y, and Z, with 
the coefficients in the above table. The calculated field changes very 
gradually, suggesting that the major part of the field represented by the 
third degree polynomial model should be originated from the earth's core. 

If the wavelength of the field is estimated very roughly by a formula 
presented in a previous section, the third degree term is to represent the 
magnetic field wavelength of 2500 km. In terms of spherical harmonics, 
this corresponds to degree n=16. At this degree, the geomagnetic field 
observed at the surface of the earth is considered to be dominated by the 
crustal fieM. One may say that the model expressed in a third degree 
polynomial series represents mostly the core field, although it contains 
small fractions of contribution from the crustal field. 

3. Comparison of the regional model with that derived from MAGSAT data 

The regional field model expressed in a polynomial series was compared 
with MGST(4/81) model (Ueda et al.. Memo. Al-2) . Field values at the grid 
point of every 1° of longitude and latitude were synthesized from two models, 
the regional field model and MGST(4/81). The difference of the regional 
model minus MGST(4/81) is shown in Fig. Al-6 for the total intensity. The 
difference is small; its absolute value is less than 20 nT almost all over 
the Japanese islands, except for the area near the border of the surveyed 
area, where the difference becomes -40 nT. Comparison was also made for the 
vector data. Discrepancy is insignificant in the central part, less than 
20 nT, for X and Z components. Regarding Y component, however, it becomes 
-80 nT; this large discrepancy might be caused by poorer orientation of 
airborne magnetometers. Hence, it may be said that, except for Y component, 
agreement between the regional model and MGST(4/81) is satisfactory. 

Similar comparison was made with the field values computed from IGRF 
1980.0. The difference is shown in Fig. Al-7 for the total intensity. 

Although a zero line passes through the central part of the surveyed area, 
the difference becomes as large as -80 nT toward the west, which is much 
larger than in the case of MGST(4/81). 

Better approximation of the model MGST(4/81) in comparison with IGRF 
seems to stem from the difference in the number of spherical harmonic terras 
involved (up to the degree and order of 13 in MGST(4/81), and 10 in IGRF 
1980.0). As was discussed earlier, the regional model contains fields with 
wavelengths of 2500 km, which roughly corrresponds to degree 16 in terms of the 
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spherical harmonic series. The IGRF 1980.0 model seems insufficient to 
approximate the field of this wavelength. It is perhaps due to the same 
reason that the past IGRF models (IGRF 1965.0 and 1975.0) have given rather 
poor approximation to the geomagnetic field over Japan and its vicinity. 

While the past IGRF models have always given 100 to 300 nT higher 
values for the total intensity than observation over the area of the Japan 
islands, the MGST(4/81) model has been confirmed to approximate the observed 
field better. The field values along the flight tracks shown in Fig. Al-4 
were synthesized from MGST(4/81), and then subtracted from the observed 
aeromagnetic data, of which running means over 18 km intervals were taken. 
Residuals computed in this way are shown in Fig Al-8, where positive and 
negative residuals are depicted in the opposite direction from the track 
line. It is seen that positive and negative residuals are distributed with 
almost equal frequency, implying that the field value synthesized from 
MGST(4/81) is close to the mean field of the observed data. 

A-2. Local Magnetic Anomalies and Their Origin 

The objective is to extract local magnetic anomalies and to investigate 
their origi.'. . There are abundant data so far accumulated by the surface 
surveys, not only on land, but also on the sea surface. Along with the 
MAGSAT data, they were originally planned to be integrated into a single 
data set to generate the local magnetic anomalies, ror this purpose it 
is necessary to reduce all the data to those at a common altitude. A 
downward continuation program by i'se of a Fourier analysis technique has 
been completed, but not yet used practically for extraction of the short 
wavelength anomalies. Only the method will be described here. Accordingly, 
MAGSAT date: were analysed independently from the surface data, and reported 
separately . 

Two approaches were made to extract the local anomalies. The MGST 
(4/81) model was subtracted from the observed values. and residuals were 
computed (Nakatsuka and Ono, Memo. A2-1). The second approach is to apply 
a polynomial series approximation to a limited segment cf the MAGSAT data 
and to compute the residuals (Nakagawa and Yukutake, Memo. A2-2). Both 
analyses have revealed interesting features in the crustal anomaly, but 
many of the detailed structures still remain unresolved. A three-dimensional 
Fourier analysis technique was explored with local Cartesian coordinates 
(Nakagawa and Yukutake, Memo. A2-2). A computer program for the downward 
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continuation has been tested with lower order harmonics up to 3 , only 
giving long wavelength features. 

Total intensity data in the oceanic area surrounding the Japanese 
islands were compiled and analysed (Fujita et al., Memo. A2-3) . Data 
acquisition covers such a long period from 1961 to 1979, that correction of 
the secular variation is needed to compare the data of different epochs. 
Definitive International Geomagnetic Reference Field (DGRF) models were 
used for the correction. The analysis is still under way. The results 
so far obtained are consistent with those of MAGSAT analysis. The surface 
data reinforce new findings made only vaguely by the analysis of MAGSAT 
data . 

In' estigation of the origin of the magnetic anomalies is the most 
retarded part of the project. Equivalent source procedure has been employed 
to infer the distribution of magnetization in the crust (Yanagisawa et al., 
1982) . There are several factors conceivable for causing the observed 
magnetic anomalies. However, the preliminary study seems to suggest that 
either the varying thickness of the crust or its thermal structure plays 
an important role in generating long wavelength anomalies. 


1. Local magnetic anomalies derived from MAGSAT 
1-1 . Anomalies as the residual field from MGST(4/81) 


Anomaly maps of total intensity and vector components were drawn for 
the residuals of the observed field values from those computed by MGST(4/81) 
model (Nakatsuka and Ono, Memo. A2-1). 

(a) Data selection: Only half orbit (ascending or descending) paths 

passing through the area in the vicinity of Japan (15°-60°N, 115°-160°E) 
were picked-up. In order to eliminate magnetic disturbance effect, the 
data for the period of K^i2 + were excluded. Attitude quality words 
(refer to NASA Technical Memorandum 82160) were checked and data points or 
poor quality in attitude (^ 4000) were excluded. 

(b) Data processing : Every orbital series of data points was tested if 

it had sufficiently uniform distribution ir. dipole latitudes, and non- 
uniform ones were dropped. After the subtraction of the reference field 

of MCST(4/81) , the magnetic field data still included external (magnetospheric 
and ionospheric) disturbance fields. A trigonometric function form 
[a • sin <{> + b • cos $ + c] of dipole latitude ( 41 ) was fit within the dipole 
latitude range of 15°-55°N, and the residual between them was assumed to 
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be the crustal anomaly field. This operation was applied to each component 
and total intensity. A statistical process is then required to draw a 
two-dimensional anomaly map. We take a local Cartesian coordinate system 
around Japan, the origin of which is located at the point of 37.5°N (geo- 
centric), 137. 5°E, and at a geocentric distance of 6800km. In order to find 
mesh point value a 3-dimensional linear regression formula was applied to 
all data within 160 km in horizontal separations f"om the mesh point to be 
reduced, and the mesh point value at the mean elevation of these points 
was determined. This regression analysis was conducted at every AO km 
mesh point, to get 2-dimensional mesh distribution of the magnetic field. 

(c) Results: The data processing described above was applied to four kinds 

of data subsets, namely A(ascending dusk-side paths only, 7587 points), D 
(descending dawn side paths only, 10A88 points), H (high-altitude paths of 
>6800 km geocentric distance, 9745 points) and L (Low-altitude paths of 
<6500 km geocentric distance, 8330 points), and the set of all available 
data (namely B, total of 18,075 data points). Fig. A2-1 shows the total 
intensity (F) anomalies around Japan extracted from each data subset A, D, 

H and L. Magnetic anomaly maps F,X,Y and Z extracted from set B of all 
available data are shown in Fig. A2-2. 

(d) Discussion 

We see in Fig. A2-1 that, in respect to the large scale anomaly, similar 
patterns are obtained irrespective of the data subset. This implies that 
the reproducibility of anomaly patterns of this scale is good. Also, the 
amplitude of anomalous field ac lower altitude (d) is larger than at higher 
altitude (c). This result is consistent with the fact that the average 
elevational difference between these two data subsets is about 100 km, if 
we consider the crustal anomaly decreases its amplitude in proportion to 
the inverse square to cube of the distance. 

On the other hand, there are some minor disturbances with shorter 
wavelength in Fig. A2-1, but they are not coincident with one another in the 
results from each at the data subsets. As the data processing step effectively 
includes a low-pass f itering operation, this process may be inadequate to 
reveal smaller- scale nomalies of shorter wavelength of a few hundred kilo- 
meters or less and also there may be the problem in the retrieval of the 
external field. Anyway, we cannot distinguish reliable magnetic anomalies. 

The anomaly maps for the three components have also been made; they 
are more noisy than the total intensity . However, if all the available data 
are analysed together, the noisy character is very much reduced as seen in 
Fig. A2-2(b)-(d). Tt is concluded that the anomaly patterns with wavelengths 



longer than about 500 km can be sufficiently reproduced by the present 
analysis . 

1-° . Anomalies computed by the polynomial series 

It has been confirmed that polynomial approximation for a limited area 

is useful to separate the crustal field from the core field (Nakagawa and 

Yukutake, Memo. A1--1) . An anomaly map was drawn for the area of 18 # -58°N 

and 120°-160°E, by means of the data on ail the satellite paths passing 

through the area of 8°-68°N and 110°-170°E for the magnetically quiet 

period K <2 . Since it is knoim that the core field becomes comparable 
p= o 

with the crustal field at degrees from 5 to 7, the 5th degree polynomials 
of latitude were used to subtract the core field, and the residuals were 
assumed to represent the magnetic anomaly of crustal origin. Anomaly 
charts were drawn at an altitude of 430km, chat is an approximate mean of 
the satellite paths. A chart of AZ anomaly is shown in Fig. A2-3. Besides 
short wavelength anomalies, some of which are still seriously contaminated 
by noisy data contained in the original data set, an east-west trend of 
long wavelength is noticeable. Two causes are conceivable for this trend. 
One is the nature of the filtering process. Although the polynomial fitting 
process filters out the long wavelength anomalies along the orbital pr.th, 
i.e. the anomalies in the north-south direction, those with a strike in 
the east-west direction perpendicular to the path remain unfiltered. 
Secondly, there is a possibility that the core field is not completely 
removed by the 5th degree polynomials, and that the core field contained in 
the 6th degree term dominates over the crustal field expressed by the rest 
of the polynomial series. 

AX, AY and AZ residuals were computed into gridded values at 1.25° 
intervals in longitude and latitude without altitude correction, and 
expressed in a double Fourier series of longitude and latitude. In order 
tc filter out the east-west trend appearing in Fig. A2-3, the first order 
harmonics were excluded from the series for synthesized. The results are 
shown in Figs. A2-4(a)-(c), where the series was synthesized up to the 
order of 31. It is clear in Fig. A2-4(a) that the east-west trend is 
almost completely removed. Instead, a weak tendency appears running in the 
northeast-southwest direction. It is not known at the present whether this 
is partly of artificial origin or not. Some are real, such as a positive 
belt along the Kuril islands from Kamchatka to Hokkaidfi. Although there 
are many anomalies of short wavelength, anomalies smaller than 500 km 
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are supposedly artifacts created by erroneous data. 

1-3 . Results derived from MAGSAT data 

A positive anomaly of total intensity and the vertical compcr.cnt seems 
to exist along the trench. An anomaly along the Kuril islands is particu- 
larly intense (Figs. A2-2(a) and (d), Fig. A2-4(a)), which corresponds to 
a narrow belt of positive anomaly observed at the surface ru. .ling along the 
Kuril Islands on the land side of the trench. On the other hand, another 
positive belt that is alsc remarkable at the surface in the north-south 
direction along the east coast of Northern HonsyQ is not clearly recognizable 
at the satellite altitude. This may suggest that the positive belt of 
Northern HonsyQ is of shallow origin in the crust, decaying rapidly with 
height, while that of the Kuril Islands has deeper roots. 

It is also interesting to point out that, over the Izu-Mariana trench, 
a positive area extends north-south (Figs. A2-2(a) and (d) , Fig. A2-4(a)). 
This agrees well with the result of analyses of the marine data as will be 
d.'scussed later, where the positive anomaly runs north-south in the vest cf 
the trench. These results seem to imply that it is a general feature of the 
magnetic anomaly in the subduct ion zone tlv _ the positive anomaly runs 
parallel to the trench on its island side If this is true, it would 
become an important key to disclose an evo. iou process of the island-arc 
and the trench structure. 

It has been believed from the surface survey that the Japan Sea is 
covered with a negative anomaly of AF and AZ. However, Figs. A2-2(a), (u) 
and Fig. A2-4(a) indicate that only the northeastern part of the Japan Sea 
is negative, whereas the southwestern half is rather positive. This may 
be related to the fact that the oceanic basin extends in the northeastern 
part of the Japan Sea, while in the southwestern part there are submarine 
rises such as the wide Yamato-tai rise. A similar feature is seen in 
Okhotsk Sea, where negative anomaly is more conspicuous. This seems to 
suggest an important nature of the back-arc basin. 

Another feature to be noted is a negative anomaly of AZ that covers the 
Korean peninsula. This forms a pair with a positive anomaly in the south. 

In the AY anomaly chart of Fig. A2-4(c) ; we may notice a quadrant structure 
at the place of the pair of AZ anomalies, where AY is positive in the north- 
east and southvest quadrants, and negative in the northwest and southeast 
quadrants. These a typical patterns of the anomalies caused by dipole 
beneath the ground with its magnetization in the direction of the present 
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magnetic field. AX anomaly is also consistent with this model. This 
suggests that a highly localized magnetic material with extremely intense 
magnetizatioi is buried, perhaps at a shallow depth near the southwest tip 
o f the peninsula. 

2 . Three-dimensional Fourier analysis of MAGSAT data 

The treatment in the previous section ignored changes in the field 
values caused by the change in the satellite altitude from 350 km to 500 km 
even on a single path. This is one of the causes generating noises that 
conceal the detailed structure of the anomalies. The altitude correction 
is an important step to refine the anomaly map. 

In section A2-1-2, the field values were expanded into a Fourier 
series of longitude and latitude. This is a practical method when the data 
are given by the geodetic coordinates. However, this cannot be used for 
the altitude correction, because the Fourier series in longitude and 
latitude do not satisfy the Laplace equation in a rigorous sense. 

in this section, we employed a local Cartesian coordinate system with 
its origin at 35°N, 140°E. The z-axis is taken vertically downwards, and 
the x- and y-axis toward the north and the east respectively on the plane 
tangent to the earth. The coordinates of the satellite position were trans- 
formed into the local coordinates, and the observed field was decomposed into 
X, Y and Z components or the Cartesian coordinates. The dc imposed values 
were then subjected to a three-dimensional Fourier analysis (Nakagawa and 
Yukutake, Memo. A2-2). 

Among the data on the paths over the area of 8°-68°N and 110°-170°E, 
only those with high accuracy of satellite attitude control were selected, 
although this process resulted in reducing the data set to about a quarter 
of the original size. The data thus selected were subjected to the polynomial 
fit. Since the core field is likely to remain in the residuals even after 
the trend was removed by the 5th degree series, we employed the 7th degree 
polvnomials to compute the residual anomaly field. Of the residuals obtained 
in this way, those inside the space -L/2<x<L/2, -L/2<y<L/2 were Fourier 
analysed, taking the maximum order of harmonics as 3. 

From the Fourier series determined in this way, it has become possible 
to compute the anomaly field at any point in the three-dimensional space. 

The anomalies at the altitude of 300 km were svnthesized and are shown in Fig. 
A2-5. Since only the lower harmonics have been used for the analysis, the 
anomaly charts delineate only the long wavelength featu • . Upward and 
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downward continuation was made and the magnetic anomalies at the altitude 
of 100 km and 500 km were computed, of which the anomalies of the vertical 

component are shown in Figs. A2-6(a)and (b) . 

Even in this crude approximation with only the lower harmonics, one may 
see the magnetic features revealed in the previous section. AZ anomaly in 
Fig. A2-5, for example, shows that the northeastern part of the Japan Sea 
is weakly negative, while the southwestern part is positive. A positive 
anomaly is noticeable from the Kuril Islands to Hokkaido. A positive anomaly 
that runs along the Izu-Mariana trench is not well recognized in this figure. 
In order to see this relatively short wavelength anomaly, it is necessary to 
include higher harmonics in the series. On the other hand, a negative 
anomaly covering the Korean peninsula is very remarkable in the present 
case. Another feature to be mentioned is, though weak, a trend in the 
direction of ENE tc WSW of alternating positive ana negative anomalies. 

Further investigation is needed to confirm whether this trend is real or not. 

3 . Local anomalies de ^ /ed from maritime data 

A large amount of daia has already been accumulated from the surface 
surveys, which can be compared with the MAGSAT data. Unfortunately, 
however, these surface data are not well compiled nor analysed in a form 
ready to compare with the satellite data except for a very few works (Nomura, 
1979). Among the data accumulated in the Oceanographic Data Center of 
Hydrographic Department, Martitime Safety Agency, total intensity data 
obtained by ship-borne proton precession magnetometers covering the region 
of 120°E-160°E and 15°N-50°N were analyzed (Fujita et al., Memo. A2-3). 

The total number of data amounts to about 600,000, which had been acquired 
during the period from 1961 to 1979. 

Since the observation covers such a long period a certain kind of 
correction for the geomagnetic secular variation is necessary to reduce the 
data to a common epoch for computing the magnetic anomalies. For the 
correction, we used the Definitive International Reference Field (DGRF) for 
1965.0, 1970.0 and 1975.0, together with the International Geomgnetic 
Reference Field (IGRF) for 1980.0. In the first place, comparison was 
made of these reference field models with the observed fields at Japanese 
magnetic observatories, Memambetsu (MMB), Kakioka (KAK) , and Kanoya (KNY), 
as shown in Fig. A2-7. The reference field values were linearly interpolated 

between 1965-70, 1970-75, and 1975-80. It is obvious from Fig. A2-7 that 
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the total intensity of the reference field is larger than that observed in 
Japan by about 100 nT, and the differences are almost constant for the two 
periods 1965-70 and 1975-80, but from 1970 to 1975 the differences increase 
abruptly by more than 50 nT at Memambetsu and at Kakioka. Therefore, the 
adopted reference field may be taken to represent the secular variation in 
total intensity in Japan with an accuracy of about 50 nT. Since the ampli- 
tude of magnetic anomalies in total intensity around Japan in about 300 nT, 
the magnetic anomaly chart based on these reference fields may not be 
seriously influenced by the correction of the secular variation. 

The anomaly was obtained by taking the difference between the observed 
values and those synthesized from the reference models. At each observa- 
tion point, the difference was computed, i.e. observed total intensity minus 
reference total intensity, and averaged every 0.25° xO. 25° area. Using all 
the data from 19b5 to 1979, Fig. A2-8 was depicted on the basis of these 
av< raged anomalies. 

For the two stable periods of the secular variation, i.e. the periods 

1965- 70 and 1975-80, the anomalies were computed separately and compared. 
Although there is a tendency that the anomalies computed from the data in 

1966- 70 have larger values than those in 1976-79, the general patterns of 
the anomalies are similar in each period. In Fig. A2-8, several features 
are remakable. (1) There are positive anomaly belts running parallel with 
trenches on their continental side, such as the anomaly along the Kuril 
trench, and tha* - associated with the Japan trench, and the anomaly parallel 
with the Izu-Mariana trench. As has already been described, the anomalies 
along the Kuril trench and the Izu-Mariana trench were visible at the satel- 
lite altitude, whereas the anomaly running north-south in parallel with the 
Japan trench was not recognizable in the MAGSAT data. (2) The southwestern 
part of the Japan Sea is covered with a positive anomaly, and an area of 
negative anomaly spreads to the north. This feature was confirmed by the 
MAGSAT date. However, the marine data suggests existence of a positive area 
again further to the north. (3) A strong positive anomaly is clearly seen 
near 40°N and 155°K. This has not been confirmed by the satellite data. 

On the other hand, the positive anomaly region in the Pacific Ocean 
marked by "?" is not clearly established. In the Pacific Ocean It is 
difficult to draw a clear anomaly pattern, due either to complex super- 
position of small scale magnetic anomalies or to sporadic distribution of 
the regions where data are missing. The number of data used for calculation 
of averaged anomaly in each unit area amounts to several hundreds near 
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Japan but it is only several in the Pacific Ocean far off Japan. 

4 . P ossible causes of the m agnetic anomalies 

Applying an equivalent source procedure, distribution of magnetization 
in the earth's crust was determined to generate the observed total intensity 
anomalies at the MAGSAT altitude (Yanagisawa et al., 1982). In comparison 
with the other geophysical data, possible causes of the anomalies were 
investigated. Lateral variation in the thickness of the crustal layer is 
one possible cause. The other is difference in geothermal structure 
revealed by heat flow measurements. These are the most likely sources for 
the magnetic anomalies. 

By placing magnetic dipoles at the surface of the earth with their 
axes parallel to the main field calculated from the MGST(4/81) model, it was 
attempted to interpret the anomaly of total intensity at a satellite 
altitude shown in Fig. A2-9, which was derived by the method described in 
section A2-1-1 by subtracting the MGST(4/81) field from the observation. 

The dipoles were placed at every 3°x3° latitude-longitude mesh points. The 
magnetic moment of each dipole was determined by means of an equivalent 
source procedure proposed by Mayhew (1979) so as to fit the observed anomaly. 
The magnetic moment thus obtained was then converted to a magnetization of 
the earth's surface layer on an assumption that the layer has a uniform 
thickness of 30 km. Fig. A2-10 shows the distribution of magnetization 
thus obtained. 

A zone of high magnetization runs almost north-south near the center 
of the Japan Sea. Most of this high magnetization roughly coincides with 
the location of the Yamato-tai, a sea platform with a shallow depth of less 
than 2,000 m. Except for this high, the magnetization in the Japan Sea area 
is lower than that on the Japan Island by 1 x 10 to 2 x 10 A/m. The weak 
magnetization is most pronounced in the north-eastern part of the Japan Sea, 
i.e., most of the Japan Basin and the Yamato Basin where the depth of the 
bottom is typically 3,000 to 3,500 m. 

Fig. A2-10 was obtained by assuming a layer of 30 km thickness for the 
source of magnetic anomalies. The anomalies in this figure can therefore 
be interpreted as representing (1) the lateral variation of susceptibility 
in the magnetized layer. There is however, another possibility for explain- 
ing these anomalies, i.e. (2) lateral variation of the thickness of the 
magnetized layer. If the thickness of the layer is assumed constant, 
variation of its thickness causes the magnetic anomalies. From comparison 
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with other geophysical data in this area, the second interpretation appears 
more plausible. 

It was tentatively assumed that the crustal meterial had a uniform 
-3 -4 

susceptibility of 4x10 in SI unit (®3.2xlO emu/cc). Then the crustal 
material is magnetized by 1.6x10 ^A/m (1.6x10 ^emu/cc) per unit volume 
under the geomagnetic field of 40 A/m (=0.5 Gauss). The magnetization 
contrast in Fig. A2-10, where the thickness of the layer is assumed to be 
a uniform 30 km, was converted to the differences of the thickness of the 
magnetized layer with the uniform magnetization, 1.6x10 ^A/m. The contrast 
of 1.3x10 ^A/m between the Japan Sea and the Japan Island in Fig. A2-10 
corresponds to the 24 km difference in the thickness. This value certainly 
depends on the assumed susceptibility, but suggests that the magnetized layer 
is substantially thinner in most parts of the Japan Sea. 

From the seismic data, crustal thickness is about 30 km beneath the 
Japan Island and about 10 km beneath the Japan Sea (Yoshii, 1979). If the 
Curie point isotherm is deeper than the crust/mantle boundary, the crust 
would be magnetized, but the mantle would not be magnetized at all 
(Wasilewski et al. , 1979). If this is the case in this area, the thin 
magnetized layer in the Japan Sea may be explained by the thin crust. 

On the other hand, the heat flow through the sea floor of the Japan 
2 

Sea is about 100 mlJ/m , and about twice the value observed on the Japan 
Island (Yoshii, 1979). The Japan Sea area is characterized by low surface 
wave velocities, which means a hot upper mantle in this area (Abe and 
Kanamori, 1970). Also the large negative residual gravity anomalies indicate 
that the upper mantle is light and hot (Yoshii, 1972). These facts suggest 
that the Curie point isotherm is shallow beneath the Japan Sea and only the 
shallower part of the crust could be magnetized. 

From the he3t flow values and geothermal gradients measured at the 
surface, the depth of the 500°C isotherm is estimated to be about 30 km 
beneath the central part of Japan (Watanabe, 1968). If this can be regarded 
as the thickness of the magnetized layer in central Japan, the 24 km differ- 
ence of the thickness leads to a 6 km magnetized layer beneath the Japan 
Sea. This is smaller than the Japan Sea crustal thickness (about 10 km). 
Therefore, this appears to suggest that the lower crust is not magnetized 
because of high temperature. 

Study on the origin of the magnetic anomaly is still at its initial 
stage. Although two highly possible causes were investigated, there is no 
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decisive evidence to determine whether the anomaly is created by the lateral 
variation of the thickness of the crust or by that of the depth of the Curie 
point isotherm. Further investigation is needed. 

A- 3 . Crustal Structure in the Antarctic 

The ultimate goal is to clarify the crustal structure in an area near 
the Japanese Antarctic base stations, Syowa Station at 69°00'S, 39°35'E 
and Mizuho Station at 70°42'S, 49°20'E, incorporating the geomagnetic data 
with gravity and seismic data. The aeromagnetic survey was the central 
program related to the present investigation. Under severe conditions the 
aeromagnetic surveys were conducted and the data are now being analysed 
(Kaminuma and Shibuya, Memo. A3-1) . 

With a proton precession magnetometer, the aeromagnetic survey was 
intensively carried out around Syowa Station (Fig. A3-1) in 1980 and 1981. The 
survey was carried out along the 20 flight lines over the eastern and central 
part of Liitzow-Holm Bay, 3 hours each, with the flight elevation of 1500 feet 
over the ground surface in 1980. The survey was also carried out in the 
Yamato Mountain area, which was located 200-300 km south from Syowa Station, 
and in the area between Syowa and Mizuho Stations. 

Fig. A3-1 shows the course of the aeromagnetic survey at 5000 feet which 
was carried out in 1381 and the total flight distance was 1685 km. Surveys 
at 1500 feet and 3000 feet were also carried out and the total flight distance 
was 1120 km and 1325 km, respectively. The total intensity of the geomagnetic 
field for three different flight elevations of 3000, 5000 and 8000 feet is 
given in Fig. A3-2. These surveys were carried out along the longitude of 
39°25'E from 68°50'S to 69°40'S in latitude. As the surveys were done within 
a few days, during a magnetically quiet period when K-index was 0-1, no 
correction has been made for the intensity values in Fig. A3-2. 

Secular variation of the geomagnetic field at Syowa Station, is given in 
Fig. A3-3. The open circles indicate total intensity. The solid line 
indicates second-order least square regression of the observed values, where 
standard deviation is 37 nT. The shorter broken line gives secular variation 
calculated by the model IGRF(1980). The longer broken line gives similar 
secular variation by 11GST(4/81). 

Using the data obtained above, an isomagnetic chart of the total 
intensity is being made, and the analyses are in progress. 
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B. ELECTRIC CURRENTS AND HYDRQMACNETIC WAVES IN THE IONO S PHERE AND 
MAGNETOSPHERE 

The objective is to investigate ionospheric and magnetospheric 
contributions to geomagnetic variations, field-aligned currents, geo- 
magnetic pulsations and hydromagnetic waves by analysis of MAGSAT data, 
and if necessary also through comparison with ground magnetic variations. 


B-1 . Ionospheric and Magnetospheric Contributions to Geomagnetic Variations 

The MAGSAT data were proved to be very useful in studying the ionospheric 
and magnetospheric contributions to geomagnetic field variations from various 
standpoints. The main results obtained so far are summarized below. 

1 . External field correction for satellite magnetic data 

A new conventional method was introduced to learn the features of 
external and internal contributions to the data obtained by MAGSAT 
(Yanagisawa and Fukushima, Memo. B 1 — 1 ) . Let AB denote the difference 
between the magnetic field B observed by MAGSAT and that given by the 
MGST(4/81 ) model. AX mag and AZ mag of the residual field in geomagnetic 
coordinates were subjected to a harmonic analysis with Legendre functions 
in the latitude range of *55°, by means of the equation 

W = a n | i {E n (^) n + I n ( 7 ) n+ '} P n (cos0) (1) 

where W is the magnetic potential for AB. The analysis was carried out 
for the dawn and dusk sides separately. The calculated values of external 
and internal coefficients (E n and l n ) of the magnetic potential W are shown 
for four classes of Dst values in Fig. B1-2 with the total number of data 
used at the right shoulder of each diagram. 

It is clear in Fig. B1-2 that E^ increases with the absolute value 
of Dst, whereas E 2 -E & are negligibly small in comparison with E^ in all 
cases. This means that the magnetic field of origin external to the MAGSAT 
level is well approximated by a simple southward magnetic field along the 
geomagnetic dipole axis. A noticeable dawn/dusk asymmetry of this southward 
field is discussed in the next section. 

The internal component of B consists of the induction current within 
the earth and the electric current flowing in the ionosphere below the 
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MAGSAT level. Ij seems to increase in approximate proportion to Ej, but 

1, -1^ do not change so much with value. We may assume here that the 
magnetic field originated from E^ and terms represents the ring current 
field and its induction effect, and the remaining fields (expecially those 
from 1,-1^) the ionospheric current as a first approximation. 

Fig. B1-3 shows the worldwide Z-anomaly map for (a) the dawn side and 
(b) the dusk side, after removing the MGST(4/81) model and the contributions 
from Ej and Ij terms. The two maps for the dawn and dusk sides differ consid- 
rably from each other because of the contamination from the ionospheric 
current. The latitudinal profiles of Z-anomaly in Fig. B 1 — 3 were averaged 
for all longitude range, and the average is shown as delta-Z in Fig. B1-4 
with standard deviations. (Fig. B1-4 confirms the peculiar latitudinal 
dependence of delta-Y found first by Maeda et al. (1982), but the discussion 
in this section is confined only to delta-Z for the dawn and dusk meridians.) 
The mean delta-Z curve in Fig. B1-4 is thought to represent the contribution 
to Z along the dawn and dusk meridians from the ionospheric current dependent 
on local time. Subtracting delta-Z of Fig. B1-4 from the Z-anomaly maps 
of Fig. B1-3, we obtain the maps shown in Fig. B 1 — 5 , in which the two maps 
for the dawn and dusk meridians are nearly the same as each other. 

In conclusion, it is shown that the magnetospheric contribution is 
we 1 1 represented by a simple southward field with a noticeable dawn/dusk 
asymmetry. The elimination of the ionospheric contributions (coming from 
below the MAGSAT level) gives a common anomaiy map in the dawn and dusk 
regions . 

2 . Dawn-dusk asymmetry of the ring current field 

The values of coefficient Ej (which represents the southward magnetic 
field of magnetospheric origin) in eq. (1) of the previous section are 
calculated for 26 days from November 2 to 27, 1979, and shown in Fig. B1-6 
(Yanagisawa and Fukushima, Memo. B1-2), together with the simultaneous 
Dst and Kp values. The E^ values generally vary in close correlation with 
Dst, but they exhibit sometimes pronounced dawn/dusk asymmetry (enhancement 
at the dusk side, and a reduction at the dawn side), especially when Dst 
and Kp values are great. This tendency is seen in Fig. B1-7, which 
demonstrates that the dawn/dusk asymmetry of E^ values is noticeable for 
Kpi3_, but almost absent for Kp^2 + . 
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The quantitative relationship between and Dst values on the dawn 
and dusk sides is shown in Fig. B1-8 with the empirical formulae, which 
demonstrate the dawn-dusk asymmetry of E^-Dst relationship but with nearly 
the same offset value of 22-23 nT for Dst=0. 

Since the ground Dst contains also the contribution from the induced 
electric current within the earth, which is a main part of 1^ in eq . (1), 
the comparison of E] + I| versus Dst was also made and is shown in Fig. 

B 1 —9 . From this figure we see that the average empirical formula for the 
dawn and dusk sides is E^ + 1^ = -1.0 Dst +■ 24. 

A comparison of the results observed by 0G0-5 (Sugiura, 1973) and 
MAGSAT was made to infer the ring current flow in the magnetosphere. OGO-5 
observed the minimum fielu depression at 2. 3-3. 6 Rg. The empirical formulae 
are 

E 1 (at OGO-5 level) = -0.83 Dst + 45 nT, 

E 1 (at MAGSAT level) = -0.74 Dst + 22 nT, 
so that E ] (MAGSAT) = 0.89 E, (OGO) - 18 nT. 

The factor 0.89 in the last equation will be attributable to the difference 
in geocentric distance of the observation level of MAGSAT (<v1 Rg) and OGO 
(2. 3-3. 6 Rg) . If we represent the ring current by a simple circular line- 
current, the factor 0.89 is obtained if the line current is assumed to 
be situated around 8 Rg. On the other hand, the difference of 1 8 nT between 
Ej (MAGSAT) and 0.89 Ej(OGO) will be reasonably attributed to an eastward 
ring current, which must exist in the inner edge of the westward ring current 
region to keep the plasma pressure of tripped charged particles in the 
magnetosphere. If the inner eastward current is also represented by a 
simple line-current at 2 Rg, the current intensity is 3 , 10'* A, which 
is not an unreasonable value (Yanagisawa and Fukushima, Memo. Bl-3) . 

3 . Electric current thr o ugh the plane encircled by the MAGSAT orbit 

It is possible to calculate the total amount of electric current 
passing through the plane enclosed by the MAGSAT orbit by means of Ampere's 
theorem (or integrated form of a Maxwell equation) shown in Fig. B-10, where 
ds is a differential segment along the satellite orbit, and B t is the 
magnetic field component tangential to the satellite orbit. The electric 
current will be reasonably thought to flow in a way shown in the figure, 
i.e. as field-aligned current in the magnetosphere, and as horizontal 
current in the ionosphere. 
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Attention must be paid in this analysis to the fact that MAGSAT 
requires 94 minutes to complete the magnetic field measurement over a 
complete circular path, and the earth rotates as much as 23.5° eastward 
under the MAGSAT orbit. This condition will result in a spurious effect 
on the calculated current under the MAGSAT orbit, unless the effect of 
the earth's rotation is corrected. 

We dealt with two kinds of Recalculations , i.e. (a) MAGSAT orbit is 

fixed in the coordinate system, the earth is rotating, and (b) the earth 

is fixed, MAGSAT orbit plane is shifting. The calculated total current J 

is denoted here as (Suzuki et al.. Memo. Bl-4) 

J(I) with MAGSAT data, and J(II) with MGST(4/81) model in case (a), 

J(III) with MAGSAT data, and J(IV) with MGST(4/81) model in case (b) . 

In the calculation in (b), the integration must be completed at the cross- 

point of the two successive loci of the satellite (cross-point between 

Tj and T. ? in Fig. Bl-11), and B^ must be taken tangential to the loci of 

MAGSAT orbits projected onto the earth (not parallel to the MAGSAT path 

in space). In a previous paper by Suzuki and Fukushima (1982), B^ was 

teken always tangential to the MAGSAT orbit in space, so that the results 

included a spurious effect due to the earth's rotation. As to J(IV), the 

integration is for a static magnetic field derivable from a magnetic 

potential, so that J(IV) must vanish theoretically. The actual result 

showed that the total current is less than 3*10^ A on a quiet day, and 

such an amount can be produced from a very small error in the absolute 

magnetic field measurement and/or in the altitude determination of MAGSAT. 

For example, the accumulation of 1 nT error in the absolute value of B 

4 c 

results in a total current of 4*10 A. An altitude error of only 150 m 
will produce a total electric current of 10 A. 

Suzuki and Fukushima (1982) published the preliminary results of their 
analysis. They compared the B^-integrations with the actual MAGSAT data 
and with the MGST(4/81) model field, as shown in Fig. Bl-12. The both 
integrals showed the total current of 4-10 6 A in magnitude with predominant 
semidiurnal variation, and the difference of these two calculations was 
assumed to indicate the electric current in space below the MAGSAT level. 
However, the authors would like to revise their conclusion in this report 
because of criticism of the method of calculation described above. 

The J(I)-J(II) calculation for a quiet day, that of November 5, 1979, in 
Fig. Bl-13 (a) in nearly the same as that of the right diagram of Fig. Bl-12 
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(which is a reproduction from the Suzuki-Fukushima paper (1982), but with a 

slight correction on the previous work). This result contains a spurious 

effect of the earth's rotation. The apparent semi-diurnal variations of 

the sunward or anti-sunward electric current in space below the MAGSAT 

2 

level is now understood mainly as an effect from g^-term of the earth's 
magnetic field (a paper for this interpretation is now in preparation). 

The reduction of semi-diurnal amplitude in the right diagram, with actual 
MAGSAT data in Fig. Bl-12, in comparison with the middle diagram, with the 
MGST(4/81) model, will be due to an induction current in the ionosphere 
below the MAGSAT level, which flows to reduce the effect of the earth's 
rotation in a conducting medium. If the ionosphere were a perfect conductor, 
the UT dependence of the curve in the left diagram of Fig. Bl-12 should 
vanish . 

On the other hand, in the results o f the J(III) calculation in Fig. Bl-13 
(b) on the same day of November 5, 1979, the total intensity of space current 
below the MAGSAT level is much smaller, i.e. less than 10^ A. We arc- 
continuing the calculation of J(ITI) for many other days during the period 
of the MAGSAT flight in order to find or check the dependence of sunward or 
anti-sunward space current below the MAGSAT level on (i) universal time, 

(ii) geomagnetic activity, (iii) solar wind and interplanetary magnetic 
field conditions, and (iv) season. At this moment, we can show only one 
example, that is on a disturbed day, November 13-14, 1979. As is shown in 
Fig. Bl-14, the space current under the MAGSAT orbit is anti-sunward in the 
disturbed condition with its magnitude far exceeding 10^ A in approximate 
proportion to the geomagnetic AE-index. 

It was described in section B-l-2 that the dawn/dusk asymmetry of 
Dst ir pronounced when the AE-index is large. Hence we may conclude that the 
increase in the AE-index is accompanied by both the dawn/dusk Dst-Asymmetry 
and the intensification of anti-sunward space current below the MAGSAT 
level. For the disturbed condition, a schematic picture of Fig. Bl-15 will 
emerge, which includes a partial ring current connected with the field- 
aligned currents in the magnetosphere and the Pedersen current in the 
ionosphere. Such a current system was proposed earlier by Fukushima and 
Kamide (1973). Although the Hall current is more important for the auroral 
electrojet, it is not shown in the schematic diagram of Fig. Bl-15, because 
the Hall current does not contribute at all to the net space current below 
the MAGSAT level because of its divergence-free flow within the ionosphere. 

It must be emphasized here that the anti-sunward Pedersen current will 
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contribute to the electrojot along the auroral oval, although the auroral 
electrojet is mainly due to the Hall current in the ionosphere. 

4 . Geomagnetic perturbations at low latitudes due to electric currents 
in the ionosphere and magnetosphere (Maeda et al. , Memo. Bl-5) 

Fig. B 1 — 1 6 shows the average perturbations (AH, AD, AZ and AF; observed 
value vinus the MGST(4/81) model) at low latitudes between 30°N and 
30°S. The difference in the observed results at dawn and dusk meridians 
is noteworthy, especially for AD, which was discussed in detail in a paper 
by Maeda et al. (1982). (The latitudinal profile in Fig. B 1 — 1 6 is similar 
to that shown in Fig. B1-4, but it must be remarked that the abscissa and 
ordinate of the two figures are drawn in opposite sense to each other, and 
the period of data used is also different. Furthermore, Fig. B1-4 eliminates 
the contribution from the ring current effect.) 

The average latitudinal profile of AD over 270°-320° meridians is 
shown in Fig. B 1 — 1 7 , where the contours of equal AD values are indicated. 

We notice here that the peculiar AD at dusk exists at altitudes lower than 
about 450 km centered at 8°N and 8”S in dip latitudes. There is another 
region of AD centered at dip latitude 5°N, at a height of about 500 km. If the 
contour lines of equal AD are approximately regarded as stream lines of 
the electric current in the ionosphere-magnetosphere, the electric current 
in the meridian plane flows across the geomagnetic field-lines indicated 
by dotted lines in Fig. B 1 — 1 7 . 

Fig. B 1 — 1 8 shows the latitudinal profile of AD at various longitudes 
and the double amplitude (ADn + ADs) is shown to be proportional to 1/B, 
where B is the magnetic field intensity at the clip equator. The altitude 
dependence of the AD range is shown in Fig. B 1 — 1 9 ; the range increases 
with decreasing altitude, from 8 nT at 500 km to 40 nT at 300 km. 

The peculiar AD on both sides of the magnetic equator seems to 
be caused by a meridional current system in the ionosphere shown in Fig. 

B1 -20, i.a. upward at the magnetic equator. The fact that the AD range 
is proportional to 1/B (in Fig. B 1 — 1 8 ) and to the sunspot number (in Fig. 
B1-21) will support the ionospheric origin of AD anomaly. Takeda and Maeda 
(in preparation) have recently carried out a calculation of F-region dynamo 
driven by pressure gradient in the evening, and obtained such a current 
system as that of Fig. B1-20. The electric, currents flowing from or toward 
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the equator are not field-aligned but flow across the magnetic field-lines. 

The meridional currents in the dusk region in the F region have never been 
found in the past, in contrast to the detection by rockets for the noon 
equatorial E region. 

A shift of the base line in AD seen in Fig. Bl-16 (positive at dusk 
and negative at dawn) may be due to field-aligned currents from the northern 
(winter) hemisphere to the southern (summer) hemisphere in the evening 
and in the opposite sense in the morning; a numerical calculation by Takeda 
(1982) for a dynamo with NS-asvmmetry seems to support this interpretation. 

Fig. Bl-22 shows the comparison of AH at dawn and dusk observed by 
MAGSAT with the ground Dst index, and Fig. Bl-23 shows a good correlation 
between che dawn/dusk asymmetry with the AE-index. (The dawn/dusk asymmetry 
with geomagnetic activity was also shown in Fig. Bl-6 and 7.). The dawn/dusk 
asymmetry of AH and its good correlation with the AE index seem to support 
the picture of current flow in the ionosphere and magnetosphere shown in 
Fig. Bl-15, i.e. the presence of partial ring current and the contribution 
of anti-su:iward current in the ionosphere to the auroral electrojet. 

5 . Su dden commencement of magnetic storms observed by MAGSAT 

From the analysis of ground nagnetograms , Araki (1977) inferred the 
presence of tlu following three kinds of electric currents flowing in the 
ionosphere during a sudden commencement (SC) of magnetic storms, i.e. (1) 
Westward zonal current induced by a sudden compression of the magnetosphere 
associated with an eastward current at the magnetopause; (2) Twin-vortex 
current for the preliminary reverse impulse on the grourd, which will be 
attributable to a dusk-to-dawn electric field in tne polar region; and 
(3) Twin-vortex type current due to a dawn-to-c'usk electric field which 
develops after SC. 

Since MAGSAT was making an accurate three-component measurement of 
the geomagnetic field above the ionosphere, it was expected to detect these 
currents flowing in the ionosphere through the comparison of MAGSAT and 
ground data, because the current flowing in the ionosphere produces a 
magne ic field in an opposite aonse at the MAGSAT level and on the ground. 

This section describes the real existence of currents (1) and (2) reported 
by Araki et al. (Memo. Bl-6), although it was partly published in a paper 
bv Araki et al. (1982). Unfortunately, it was impossible to check current (3), 
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because the MAGSAT records suffered from high-latitude disturbances 
soon after the SC's :oncerned. 

Twenty-two SC's were observed according to the Solar Ceophysical Data 
(edite by J.V. Lincoln, published by WDC-A for STP, NOAA) during the MAGSAT 
observation period of November 2, 1979 - May 17, 1980, Araki et al. (Memo. 
Bl-6) examined 21 cases of SC's out of 22 cases. 

In low latitudes, MAGSAT detected SC mainly in H-component with a 
slightly greater amplitude (average value 1,3) in comparison with the ground 
SC data. In high latitudes, MAGSAT detected SC in both H- and D-components, 
but their amplitude is i ot always measurable (especially nt latitudes higher 
than 60°) because of the contamination form polar disturbance. The larger 
amplitude of the H-component SC at the MAGSAT altitude to that on the ground 
will be due to a westward shielding current at the time of SC. The Z- 
component shows no significant change at the time of SC. 

A detailed case study w„s made for two large SC's observed on November 
30, 1979, at a low latitude, and on March 19, 1980, at a high latitude. 

Fig. BI-24 shows the SC on November 30, 1979, observed by MAGSAT at 
.-he altitude of 550 km at 26.5°S, 30.3°W over the South Atlantic Ocean 
and on the ground (Hermanus in South Africa, Trelew in Argentina and its 
conjugate point Sam Juan in Puerto Rico). This SC showed a preliminary 
reverse impulse (PRI) in D-component, and MAGSAT recorded an eastward 
deviation with its maximum magnitude of 10 nT within 3 minutes after the 
SC. A worldwide equivalent overhead current-system for this PRI is drawn 
from the ground data in Fig. Bl-25 for the northern hemisphere, assuming 
that the current flow is symmetric with respect to the geomagnetic equator. 
At Eusebio (3.9°S, 58.4°W) in Brazil (+5.6° geomag. lat.), a ground station 
close proximity to MAGSAT at the time of SC, the overhead current res- 
ponsible for AD of PRI was northward. If the symmetry of the PRI current 
is assumed, AD of PRI at the ground subsatellite point of MAGSAT must have 
been eastwarc. Since a westward AD of PRI was observed by MAGSAT, Araki 
et al. (1982) concluded that the case of SC on November 30, 1979, was the 
first experimental evidence of the global ionospheric current-system for 
PRI of SC. 

Another example of SC with PEI is shown in Fig. Bl-26, which took 
place on March 19, 1980, when MAGSAT was passing over Magadan (60.0°N, 

151. 0°E) in east Siberia. T he D-component at MAGSAT first decreased by 
43 nT in about 2 minutes and recovered to the original level. A little 
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bit later, the H-component increased by 52 nT in about 1 minute. After 
0621 UT irregular disturbances masked the SC variations. The right panel 
of Fig. HI -2b shows the time variation in the horizontal geomagnetic vector 
along the MACS AT orbit. The horizontal geomagnetic vector deviation of 
the ground PI from tin* pre-SC level is also shown in the right panel by 
dashed arrows. We see clearly that the deviations above and below the 
ionosphere are in approximate ly opposite to each other. This case also 
supports the conclusion that PR1 is caused by an ionospheric current. 


B-2. Field-AI ign ed Cu rrents 

MAGSAT provides us with the best available data in the dawn and dusk 
r *gions to studv the field-aligned current distribution, because ot its 
accurate measurement ot the three— component geomagnetic field and its flight 
at lower altitudes in comparison with other previous satellites. Insofar 
as the field-aligned currents in high latitudes arc concerned, there are 
the following thro regions of large-scale t ie ld-al igned currents, according 
to present knowledge, l.e. 

(1) Region I: average invariant latitude 72-78“ on quiet days and 68- 

75° on disturbed 'ays; current into the ionosphere on the dawn side and 
out oi the ionosphere on the dusk side; total current 3-5* H) ( ’ A without 
much influence ot magnetic activity; 

( 2 ) Regions 2: invariant latilule 65-/2" on quiet days and 62-68 on 

disturbed days, equatorial side adjacent to Region I; current out of the 
ionosphere on the dawn side and into the ionosphere on the dusk side; total 
current 2-5. 10* A dependent on magnetic activity; and 

1) Cusp Region: located poleward of an 1 adjacent to Region 1 near local 

h urn; total current 1 . 5 - 10 A; pattern is modified by B y of interplanetary 
i ignet ie f ie Id. 

In a previous paper by lijima el al. (1982), they analyzed the MAGSAT 
record during i severe magnetic storm (November I )-14, 1979). In their anuly 
sis the MGST(4/8I) model 1 ie Id was first subtracted from the actual 
MAGSAT data, and the residual was illustrated for \F, AB„ and AB a , when; 
the hitter two are Lbe components parallel and perpendicular to the main geo- 
magnetic t leld 1 Ine of the MGST (4/81) model at each observation point. Al$ A 

was further decomposed into \D (dusk-to-dawn component) and \S (sunward 

XM YM 

component), or into geomagnetic nort h-sout h (AD ) and east-west (AB ) 
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components. AB„ is very useful in discussing the contribution from electric 
currents flowing in the ionosphere of enhanced electric oncductivity during 
magnetic storms, because the MAGSAT altitude is low enough to detect the 
influence of Hall current that attenuates steeply with height. 

In this report a new result of analysis is described, which dealt 
with the field-aligned current distribution over the polar-cap region for 
an extremely quiet period when the interplanetary magnetic field was northward 
with a magnetiude as large as 5—20 nT. Such an interesting case has never 
been reported, so that the present report introduces new knowledge on the 
dependence of field-aligned current patterns on solar wind conditions. 

Fig. B2-1 shows two examples of the latitudinal profile of the residueal 
magnetic field in the northern hemisphere on Mav 10, 1980. As is usual 
with quiet days, the locations of ' wgions 1 and 2 were several degrees 
higher in comparison with disturbed days; Region 1 currents were detected 
at 75.5-79.5° invariant latitude on the dawn side and at 74.9-78.4° on 
the dusk side. The most striking feature on this quiet day is the large 
perturbations observed poleward of Region 1 with a fairly large amplitude 
of AB„ or AF. Two examples of the observed characteristics during consecutive 
paths of MAGSAT in Fig. B2-2 show the presence of upward current in the 
morning regior and downward current in the afternoon region . The intensity 
of this field-aligned current confined to the region above 80° latitude 
was 1.7 times the Region 1 current intensity, and the current direction 
is opposite to Region 1 current. The concurrent AB,, or AF observed by 
MAGSAT indicates the presence of a twin-vortex horizontal current confined 
to the polar-cap region, with an anti-sunward current over the geomagnetic 
pcle that is opposite to the ordinary S^ current (Iijima et al.. Memo. B2-1) . 

The results of analysis of ground magnetograms at 8 high-latitude 

stations also support the above conclusion. The overhead current arrows 

for the horizontal geomagnetic variation (deviation from the mean level 

between 2100-0300 MLT of May 10-11, 1980) in high latitudes in Fig. B2-3 

clearly show the existence of a reversed S^ current near the geomagnetic 

9 

pole. This is further confirmed by Fig. B2-4, which shows the simultaneous 
geomagnetic variation of AZ, namely positive in the region of a counter- 
clockwise ionospheric current, ai.u negative a the region of a clockwise 
current in Fig. B2-3. 

Although the occasional appearance of a reversed current near the 
geomagnetic pole has been advocated by Iwasaki (1971), Maezawa (1976) and 
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others from a comprehensive analysis of ground magnetograms at high-latitude 
observatories during the period of northward interplanetary magnetic field, 
the case on May 10, 1980 is the first evidence for such a current in the 
ionosphere revealed in the MAGSAT data. 

B-3. Geomagnet ic P ulsations and H vd ronagne tic Wave s 

It was planned, when the Japanese MAGSAT Team presented its proposal in 
February 1979, to study short-period hydromagnet ic waves in the earth's environ- 
mental space through spectral analyses of geomagnetic fluctuations observed by 
MAGSAT and simultaneous data on the ground. However the team members for this 
item had to fulfil other obligations during 1980-82, so the this subject re- 
mains a task for the future . We hope, however, that a number of interesting results 
will come out in the near future after comparing the MAGSAT data and other simul- 
taneous data, such as geomagnetic, VLF wave and auroral observations at various 
places in high latitudes. 
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'c) z-componcnt 


(d) total intensity F 


Fig. Al-5(a)-(d) 


Regional magnetic chart of X, Y, and F for the 
epoch of 1980. 















Fig. A2-1 Magnetic anomalies of total intensity (F) around Japan 
extracted from different data sets; i.e. (a) da a subset A 
(ascending paths only), (b) data subset D (descending paths 
only), (c) data subset H (high altitude paths), and (d) data 
subset 1. (low altitude paths). Contour interval is 1 nT. 









A2-2. Magnetic anomaly maps for respective components extracted 
from all the available MAGSAf data, for (a) *-otal intensity, F, 
(b) northward component, X, (c) eastward comprnent, Y, and (d) 
downward component, Z. Contour intervals are all 1 nT. 
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AZ ( Original ) 


A2-4(a). Magnetic anomaly AZ 
synthesized from Fourier series 
with the first harmonic excluded. 
Shaded area indicates positive 
anomalv. Contour interval is l nT 


-3. AZ residuals after removal 
of the general trend approximated 
!>y the S t h degree polynomials. 
Shaded area indicates positive 
anomaly Contour interval is 1 nT 


AX ( LowC'it) 


-4(b). Magnetic anomaly AX 
synthesized from Fourier series 
with the first harmonies excluded. 
Shaded area indicates positive 
anoirilv. Contour Interval is l nT 


Fig. A2-4(c). Magnetic anomaly AY 

synthesized from Fourier series 
with the first harmonics excluded 
Shaded area Indicates positive 
anomaly. Contour interval is 1 n 
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values from the reference field models from 1965 to 1979. Models used are DGRF 1965 

at observatories Memambetsu (MMB: 63° DGRF 1970, DGRF 1975 and IGRF 1980. 

54 ' 30"N, 144°11 ’ 35"E) , Kakiok? (KAK: 

36°13'45 ,, N, 140°11 ' 23"E) , and Kancya 
(KNY: 31°25' 14"N, 130°52 ' 56"E) . 





ORIGINAL P r .C? 15 
OF POOR QUALITY 


F-12 



Fig. A2-9. Magnetic anomaly of total intensity 
(in nT) at the average altitude of 450 km. 



longitude 

Fig. A2-10. Distribution of magnetization (in the 

unit of 10“ * A/m or 10”^ emu/cc) computed from 
the magnetic anomalies of total intensity in 
Fig. A2-9 hv use of an equivalent source 
technique . 
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Fig. A3- 1 . Geographical map around Syowa Station. Thick lines 
indicates flight lines of aeromagnetic survey conducted 
in 1981. 
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Fig. A 3-2 . Total intensity obtained by aeromagnetic surveys in 

the vicinity of Syowa Station, Antarctica, at the altitudes 
of 3000, 5000 and 8000 feet. 



Fig. A3-3. Secular variation in the total intensity of geomagnetic 
field observed at Syowa Station. 
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P-15 



W= aJ^EnCj/ 1 * In(-f) nM )P n (cose) - 


Fig. Bl-1 . Dependence on magnetic latitude of AX ma g and £Z m ag measured by 
MAGSAT for calculating the magnetic potentials of external (E n ) 
and internal (In) parts by means of the equation shown at the 
right bottom corner. 
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Fig. Bl-2. External (E n ) and internal (In) coefficient values of the mag- 
netic potential for four classes of Dst values, computed sepa- 
rately for the orbits on the dawn side and the dusk side. 
Numbers on the right shoulder of each block are those of 
orbits used for the analysis. 
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Fin. Bl-A. Latitudinal dependence of the "mean ionospheric field" on 

the dawn and dusk sides, obtained after removal of MGST(4/81) 
model field and the ring current and induction fields. The 
thl'k center curves are the average, and the thin upper and 
lower curves show the range of standard deviations. 



Fig. Bl-5 (b) . Dusk-side Z-anomaly map obtained after subtracting the mean ionospheric field 
the map of Fig. Bl-3(b). Contour interval is 2 nT. 




Dawn side data 
Dusk side data 






Fie Bl-6 Top: the coefficients E-. calculated for dawn (o) and dusk (+) paths agmst orbit number. 

Ei ‘represents the north-to-south uniform magnetic field parallel to -the geomagnetic dipole axis 
Middle: Dst values at the tim • ~f dawn equator crossing (o) and dusk equator crossing (+) of 
MAGSAT. Bottom: Kp indices at dawn (o) and dusk(+) equator crossing. 
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Bl-7. Comparison between Ej values on the dawn and dusk sides for the Kp ranges of (a) less than 

or equal to 0.^, (b) between 1- and 2+, and (r.) equal to or greater than 3-. The solid lines 
in the diagrams represent the lines of equal ' 6 \ values on the dawn and dusk sides. For high 
Kp values, the dawn/dusk asymmetry of E\ is noticeable. f 
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(a) 


Dawn side 


•ISO. DU llli.utl -70 00 90 00 

OST 


N = 2597 
R = -0.90 

El= -0 59* Dst -22 



Wj'Oo ..9.0.00 


(b) 


Dusk side 



1 r 1 1 r~ , _ . 

•I SO. 00 -110.00 -70.00 -30.00 * .$0.00 

DST 


N=2205 
R = -0.93 

El - - 0.92 xDst-23 


Fig. Bl-8. Intensity of the ring current field at the MAGSAT 
level versus Dst values on the ground. The number 
of data points (N), correlation coefficient (R) and 
the best fitted linear function are also shown. 


(a) 


(b) 



R = -0 92 

El *H =0 79 * Dst *19 


Dusk side 





0.00 - 1 10. 00* -70.00 ' - 30.00 ' ^ M S0.C , 1 

DST 


N = 2205 
R =-0.93 

E 1 * 1 1 =-1 2*Dst-29 


Fig. Bl-9. Intensity of the ring current field plus earth's 
induction and/or ionospheric field at the MAGSAT 
level versus Dst values on the ground. The number 
of data points (N), correlation coefficient (R) and 
the best fitted linear function are also shown. 
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Fiij. m-IO. Ampere's theorem Fiq. Dl-11. Two methods of Bt~ integration : 

applied to a polar orbit (a) Integral circuit in the sun-earth 

satellite. coordinate system, (b) Loci of satellite 

orbits projected onto the earth surface. 




Not At Ion 

Pm a Ini'll 

Line- Intent .it 
along th«* MA«- 

mi l*.» t 

<A7 Orbit 

1 | IIM . 

mai;| 

HAliSAl 

( i .*• 

Midi 

to 

H ? <N) 

i io«. 

MAC | 

MACSA1 

f r **• 

nm 

to 

41 1 N > 

1 |HS. 

MAC | 

MAC-SAT 

f r .hd 

H,(S> 

to 

»MS) 

1 |CS. 

NAi:j 

MAISAT 

f r vHD 

<(Sl 

t o 

Ci •;> 

I (KN. 

HI FI 

HUST(t/8l) 

f r in 

H, (Ml 

to 


III*. 

Rtri 

Mt;sT(*/Hl) 

f r .HD 

CIM) 

to 

C(N) 

1 Ins. 

mrj 

Mi;ST(A/*l) 

f f .HD 

H i ( S I 

t o 


i l< s . 

Rtn 

Mr.ST(*/«n 

f r .hd 

t\s> 




II in this fi<inri' is The same as T in 
Kiq. Dl-11, ami C denotes the cross- 
point of two consecutive orbits pro- 
tected onto the ground. N and S 
refer to the northern and southern 
hemispheres . 


Fig. nl-12. Results of the numerical line-integration of P t with the actual 

MAOSAT data (left diagram), with MOST (4/8 1 ) model reference field (middle 
diagram), and the difference between them (right diagam) . These are re- 
produced from a previous paper (Suzuki and Fukushima, 1982) . 
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Fig. B 1 — 1 3 . Total intensity of the space current through the 
plane enclosed by the MAGSAT orbit on a quiet day 
of November 5, 1979, obtained from the calculation 
of (a) J(I)-J(II), and (b) J(III). The result in 
(a) contains an effect of the earth's rotation 
under the MAGSAT orbit. The effect is eliminated 
in the calculation in (b). 
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Fig. Bl-14. Total current calculated from J(I1I) and the geomagnetic AE-index on the disturbed 
days of November 13-14, 1979. 
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Fig. Bl-15. A schematic model model of the ring current during magnetic 
disturbance' with partial ring current connected with anti-sunward 
field-aligned currents into and out of the ionosphere. This model 
explaing the dawn/dusk asymmetry of AH in low latitudes and the 
space current under the MAGSAT orbit. 


DIP LATITUDE DUSK ( 191 days, 412 km) 



AH AD AZ AF 


dip latitude DAWN ( 190 days, 419 km ) 



Fig. Bl-16 . Average profile of the residuals (AH, AD, AZ, and AF) 
from the main field at low latitudes for all days of 
the MAGSAT observation period at the dusk and dawn 
meridians, where the number of days used and the mean 
altitude of MAGSAT are shown in parentheses. 
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DIP LATITUDE (deg) 


Fig. Bl-17. Distribution of AD in the average meridian plane 
(\ = 270°-320°) , where the equal AD lines (in nT) and the 
geomagnetic field lines are shown by full and dotted lines, 
respectively. 


MAGSAT o-comp. average/ ousk-sioe 

79/12 / 02 - 12/19 20 «t 




GEOGRAPHIC LONGITUDE 


Fig. Bl-18. Longitudinal dependence of AD averaged over 18 days 

on the dusk side (top), and that of the total range (ADn+AD s ) 
together with the values of 1/B (bottom). 
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Fig. Bl-24. The time variation of the three components and the total 

force (all are the deviations from MGST(6/80) model field) observed 
by MAGSAT and the normal-run magnetograms from 3 stations near the 
MAGSAT position at the time of SC on November 30, 1979. Disconti- 
nuous jumps in MAGSAT data due to the attitude inaccuracy are 
corrected by dotted lines. 





Fig. Bl-25. Ionospheric current system for the averaged PRI in the 
northern hemisphere in the geomagnetic latitude - local time 
coordinates. Arrows indicate the direction of the N-S component 
of an assumed ionospheric current for PRI observed by MAGSAT and 
at ground stations. Positions of MAGSAT and three southern 
hemisphere stations (Trelew, Hermanus and Port aux Francais) 
are projected to the conjugate points in the northern hemisphere. 
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Fig. Bl-26 . SC event on March 19, 1980. Broken arrows in the right panel 

indicate horizontal magnetic vectors of the PI of SC observed at ground 
stations. Solid arrows show the successive variation in the horizontal 
magnetic vector along the MAGSAT orbit for the time interval shown by a 
dashed horizontal bar in the central panel. Solid horizontal bars in 
the central and left panels indicate an equal time interval. 
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Fig. B2-1. Geomagnetic vector residuals separated into AS (perturbation trans- 
verse to the main geomagnetic field and directed toward the sun) , AD (per- 
turbation transverse to the main field and in the dusk-to-davn direction) 
and ABh (perturbation parallel to the main field) , obtained from the MAG- 
SAT data for two paths over the northern hemisphere on a very quiet day ot 
May 10, 1900 (UT 1015-1045 and 2222-2252). 
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